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Monte Carlo Study of Two-Color QCD with Finite Chemical Potential 
- Status report of Wilson fermion simulation - 
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Using Wilson fermions, we study SU(2) lattice QCD with the chemical potential at (3 = 1.6. The ratio of 
fermion determinants is evaluated at each Metropolis link update step. We calculate the baryon number density, 
the Polyakov loops and the pseudoscalar and vector masses on 4 4 and 4 3 x 8 lattices. Preliminary data show 
the pseudoscalar meson becomes massive around fi = 0.4, which indicates the chiral symmetry restoration. The 
calculation is broken down when approaching to the transition region. We analyze the behavior of the fermion 
determinant and eigen value distributions of the determinant, which shows a peculiar "Shell-and-Bean" pattern 
near the transition. 



1. INTRODUCTION 

Lattice study of QCD has been expected to 
provide useful informations to understand non- 
pcrturbative aspects of quark/gluon physics. Es- 
pecially at finite temperature, it predicts the con- 
finement/deconfinment transition and is able to 
describe many features of hadrons and quark 
gluon plasma (QGP). Lattice QCD offers a sound 
base of QGP physics, which has become very im- 
portant issue of physics because of currently go- 
ing active experiments at CERN SPS and BNL 
RHIC. Sec ||. 

Contrary to the finite temperature calculation, 
the progress in lattice QCD study of the finite 
density has been rather slow. This is because of 
the well known complex action problem. Indeed 
after the first QCD dynamical quark simulation 
with the chemical potential was done for SU(2) 
color group ||], to our knowledge, no full SU(3) 
QCD calculations had been tried. A trial to put 
the phase coming from the determinant into ob- 
servables suffers from large fluctuation even at 
4 4 size lattice near the phase transition || 
Stephanov shows that the quench approximation 
is not the correct N f = limit of full QCD |§. 
We still wait for good news concerning Glasgow 



method (see a good review by Barbour j?J and 
references therein) , and finite density method Q . 

Recently the situation has been changed; Due 
to the progress in analytical investigations P,^0[ , 
we have a hope to obtain informations on con- 
cerning QCD by studying QCD-like theories such 
as SU(2) QCD, models with quarks in the adjoint 
representation and QCD at finite isospin density; 
they are expected to have less difficulties in nu- 
merical analyses. In these years, there are indeed 
high activities in Monte Carlo calculations with 
dynamical quark of these models |ll| [l^] . 

In this paper, we report our recent work on 
SU(2) QCD with Wilson fermions to study finite 
density states. 

2. ALGORITHM 

The chemical potential, /i, is introduced in the 
fermion action, ipWip, as 

W(x,x') = 5 X . X . 
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-«E{( ] - 7i)0i(aO<k )S+ J +(1 + li)U}(x')S x , iX _ 

i=l 

- K {e+^(l- 74)^)^+4 
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+e-^(l + - H )Ul(x')5 x , x _ i }(l) 

by Hasenfratz and Karsch JlTj to avoid an infin- 
ity in the energy density. As the lattice spacing 
a tends to zero, Eq.(Q) gives W(/j.) — W(0) + 
Kfia-tp^f^ + 0(a 2 ). The above formula was inde- 
pendently obtained in Ref. (r| by the following 
naive argument: In the continuum perturbation, 
the chemical potential is introduced by the sub- 
stitution p4 — » p4 — in in fermion propagators. 
In order to have this continuum limit, the lattice 
fermion propagator should have the form, 

3 

V(l -«£{(! - l*V Pia + (1 +7«)e"^ a } 

-nUl - 74 )e l ( p4 -^ a + (1 + 74 )e- l(p4 - v)a }) (2) 

Little is known about the behavior of dynami- 
cal fermion simulations when the chemical poten- 
tial is introduced. We therefore decide to employ 
an algorithm where the ratio of the determinant, 

,3, 

is evaluated explicitly at each Metropolis update 
process, U -> U + AU, where AW = W(U + 
AU) - W(U) |§-§§. An essential ingredient of 
the algorithm is the following Woodbury formula, 

(w + Awy 1 

= W' 1 - W~ X AW (I + W^AW)' 1 W- 1 (4) 

Suppose we update link variables U^xYs only 
on a subset H of whole lattice. Then AW ^ 
only on H. Woodbury formula (|J) holds even if 
the matrix space is limited on H . In this case we 
can get the ratio of the fermion determinant as 
far as U^xYs are updated inside H. We take a 
2 4 hypercube as H . When we go to the next hy- 
percube, (W~ 1 )h 's are initialized by CG method. 

3. RESULTS 

First we calculate the expectation value of the 
density, 
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Figure 1. Baryon number density divided by T 3 
as a function of \x. 



where V s is the spatial volume N x N y N z . In Fig.|l|, 
we plot < n > /T 3 , as a function of //, which 
is dimensionless. The dotted line corresponds to 
the free quark case obtained by setting (x) = 1 
and k — 1/8. When k becomes large, the density 
reaches to the free case quickly. 

The Plyakov line < L > also increases as a 
function of /i as shown in FigJ|. When the hop- 
ping parameter, re, becomes large from 0.158 to 
0.185, values of Polyakov line increase. Since 
the lattice size is small, no sharp increase of 
< L > is seen, but large values of < L > in- 
dicate quarks become free from the confinement 
force at [i > 0.4. 

The behavior of < L > suggests that we are 
near the phase transition. But the calculation 
breaks because of numerical instability at re = 
0.158, fi — 0.8 and re = 0.185, [i — 0.7 in case of 
4 4 , and we cannot go beyond. To see origins of the 
instability, we measure the ratio of the fermion 
determinant detW{U + AU)/ActW{U). Figure | 
shows the behavior of the ration as a function of 
Monte Carlo sweeps. The ratio changes the value 
around one, but suddenly it fluctuates very large 
and the calculation is broken. 

In Fig.^, we plot eigen value distributions of W 
for fj, =0.0, 0.4, 0.6 and 0.8 at re = 0.156 on 4 4 
lattice. As \i increases, the distribution of eigen 
values, Xi, becomes wide in the Real axis, and 
MinRc(A) < 0; on the other hand, A^'s scatter 
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Figure 2. Polyakov loop expectation value as a 
function of fx for k=0.158, 0.180 and 0.185. 



Figure 4. Eigen value distribution of W for 
fx=0.0, 0.4, 0.6 and 0.8 at k = 0.156 on 4 4 lat- 
tice. 
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Figure 3. The ratio of detVF as a function of 
Monte Carlo sweeps. It fluctuates around one, 
but before the simulation is crashed, big fluctua- 
tion is observed. 



sparsely around origin due to dynamical fermion 
simulation which includes detjy in the measure. 

These behaviors are expected. But the behav- 
ior of the eigenvalues around the real axis is new 
to us. As [i deviates from zero, the region near 
the real axis becomes dilute, but it seems that a 
new group appears when fj, becomes further large, 
and near the phase transition there are outer and 
inter groups like the shell and bean. This shell- 
and-bean structure occurs when the calculation 
breaks, and might be related with the phase tran- 
sition. 

We evaluate pseudo scalar and vector meson 
masses for (i =0.0, 0.2, 0.4 and 0.5. Because the 
lattice is small, we fit propagators at N t =2,3,4,5 
and 6 by one-pole fit. The result is shown in Table 



1. Although the data are still very preliminary, 
we extrapolate them to the chiral limit obtain 
Fig.^. Error bars are very large, but the pion 
mass becomes massive around /i ~ 0.4, which 
means that the chiral symmetry is restored in 
these regions. 



4. CONCLUDING REMARKS 

We present numerical study of two-color QCD 
with the chemical potential with Wilson fermions 
at j3 = 1.6. Although the lattice is very small, 
most data suggests we are reaching the confine- 
ment / deconfinement phase transition. 

We employ an algorithm which takes into ac- 
count the ratio of fermion determinant exactly, 
and has large Markov step, but we suffer from nu- 
merical instability and can not go over the phase 
transition. Near the phase transition, the dis- 
tribution of eigen values of W shows a peculiar 
"Shell-and-Bean" structure. 

Since the calculation is done at strong coupling 
region, the strange behavior of detVF might be re- 
lated with rough configurations far from the con- 
tinuum Q . We plan to continue the analysis by 
using improved gauge actions to clarify the point. 
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Table 1 

Pscudo scalar and vector masses at /i=0.0, 0.2, 
0.4 and 0.5 for k=0.158 and 0.180. 
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p >-EH 





7T 




K = 0.158 


k = 0.180 


0.0 


1.775(31) 


1.327(51) 


0.2 


1.722(09) 


1.305(37) 


0.4 


1.738(11) 


1.428(51) 


0.5 


1.649(18) 


1.346(65) 




P 


0.0 


1.882(45) 


1.509(41) 


0.2 


1.783(13) 


1.490(41) 


0.4 


1.788(36) 


1.660(49) 


0.5 


1.798(46) 


1.482(100) 
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